Energy and momentum entanglement in parametric downconversion 
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We present a simple treatment for the phenomenon of parametric downconversion considering 
the coherent scattering of one pump photon into a photon pair by a nonhnear crystaL The energy 
and momentum entanglement of the quantum state of the generated twin photons are seen as a 
consequence of the fundamental indistinguishability of the time and the position in which the photon 
pair is created inside the crystal. We also discuss some consequences of the system entanglement. 
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I. INTRODUCTION 

The generation of twin photons by a nonlinear crys- 
tal in the phenomenon of parametric downconversion has 
been extremely useful to study fundamental aspects of 
quantum mechanics and to physically implement quan- 
tum information protocols in the last decades. This 
occurs because the photon pair generated in the pro- 
cess can be highly entangled in many degrees of free- 
dom like energyji momentumyS, angular momentum^ and 
polarization.— For a pair of particles in an entangled state, 
we cannot associate an objective reality to the state of 
each of the particles, since the particle correlations can- 
not be explained using a local and realistic description of 
the system, having profound consequences in our way to 
see naturej^iS, Entanglement is also a powerful resource 
for quantum information, permitting the execution of 
quantum algorithms that are more efficient and secure 
than the classical ones.'' 

In this work we derive an expression for the quan- 
tum state of the twin photons generated by parametric 
downconversion using a simplified version of the example 
presented in our recent work,^ in which the interaction 
between light and matter is treated using the Bialinicki- 
Birula-Sipe photon wave function formalism.^' (For di- 
dactical discussions of the photon wave function formal- 
ism, see Refs. [Hil) The reader interested in a more 
formal treatment can consult Ref. H. The main idea is 
that when one photon is converted in two in a nonlinear 
crystal, there is a fundamental indistinguishability of the 
position and the instant of time in which the photon pair 
is generated. This fact forces us to coherently add all the 
probability amplitudes. The interference of the proba- 
bility amplitudes of generation at different times leads 
to the energy entanglement of the twin photons, while 
the interference of the probability amplitudes of genera- 
tion in different positions inside the crystal leads to the 
momentum entanglement of the twin photons. 

In Sec. |lT]we discuss in more detail the physical argu- 
ments of the treatment and arrive at a general expression 
for calculating the twin-photon wave function, Eq. 



In Sec. mil we compute the twin-photon wave function 
using reasonable approximations. In Se. II VI discuss some 
consequences of the system entanglement. Finally, in Sec. 
|V]we present our concluding remarks. 



II. PARAMETRIC DOWNCONVERSION AS A 
SCATTERING PHENOMENON 

If an electromagnetic wave interacts with a small trans- 
parent object that responds linearly to the electric field, 
such that the induced electric dipole moment of the ob- 
ject is proportional to the incident electric field, the resul- 
tant scattered electromagnetic field is a superposition of 
the incident field with the field generated by the oscillat- 
ing dipole moment of the object. The scattered field will 
in general have a diffraction pattern from which we can 
deduce some properties of the scatterer. If the incident 
electromagnetic field is composed by only one photon and 
we use a quantum language, we can say that there is a 
probability amplitude for the incident photon to be in- 
stantaneously absorbed and re-emitted by the object that 
must be coherently summed to the probability amplitude 
for the photon to pass directly without any interaction 
with the object in order to compute the wave function of 
the scattered photon. This instantaneous absorption and 
re-emission of a photon can be interpreted as the absorp- 
tion of the incident photon inducing the object to have 
an oscillating dipole moment and this oscillating dipole 
moment creating the scattered photon. 

For instance, the propagation of an electromagnetic 
wave through a transparent linear medium can be seen 
as a superposition of the incident wave being transmitted 
directly through the medium with the waves generated 
by the oscillating charges in the mediumJ^iii Each com- 
ponent of the resultant wave (the incident and the waves 
emitted by each scatterer) propagates at the speed of 
light in vacuum, but an interference effect makes the re- 
sultant wave to propagate with a reduced velocity. 

In the phenomenon of parametric downconversion, 
light interacts with a nonlinear medium whose nonlin- 
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ear polarization is proportional to the square of the in- 
cident electric field (or to the product of different com- 
ponents of the electric field) . A nonlinear scatterer like 
this can absorb two photons and emit one, since the in- 
duced nonlinear dipole oscillations can create one photon, 
but these oscillations occur at the expense of the absorp- 
tion of two photons from the incident field. This situa- 
tion corresponds to the phenomenon of second-harmonic 
generation. By symmetry it must also be possible that 
the nonlinear scatterer does the opposite: to emit two 
photons absorbing one, which corresponds to parametric 
downconversion. 

If we have a medium with a small nonlinear scatterer 
placed in a fixed position r', the interaction between an 
incident electromagnetic field and this nonlinear scatter- 
ing element causes a probability amplitude for one in- 
cident pump photon to be converted into two photons 
during its passage through the medium. The situation is 
illustrated in Fig. [1] where it is considered that photon 1 
is observed in the position ri and photon 2 in the position 
r2, both at time t. Of course, there is also a probabil- 
ity amplitude that no downconversion occurs (and this is 
what happens in most cases), but from now on we will 
make a post-selection to the cases in which a photon pair 
is created. We will first consider that the incident pump 
photon has a wave function ipp(r',t') in the nonlinear 
scatterer element and that the downconversion occurs at 
time t', such that the probability amplitude of generating 
the twin photons at time t' is proportional to ^/'p(r',i'). 
If we consider that both photons are born in a precise 
position r' and at a precise time t' , according to the un- 
certainty relations they must have a huge uncertainty in 
their momenta ftk^ and energies hui — hkic/ui, where 
is the wavevector, uji the angular frequency and Ui the 
refractive index for photon i (i = {1,2}) in the propaga- 
tion medium, |ki| = ki, his the Planck constant divided 
by 27r and c the speed of light in vacuum. This means 
that to compute the probability amplitude of finding the 
photons 1 and 2 in the positions ri and r2 at time t, 
we must consider the propagation of photon 1 from the 
space-time point (r',t') to the space-time point (ri,i) 
with all possible wavevectors ki and angular frequen- 
cies uJi — kic/ni, and the same for photon 2, coherently 
adding the corresponding probability amplitudes. When 
propagating from (r',t') to {ri,t), a photon in a plane- 
wave mode with wavevector k^ and angular frequency uJi 
accumulates a phase k^ • (r^ — r') — uji{t — t'), that must 
be considered in the probability amplitude. So we can 
write the probability amplitude of finding the generated 
photons in the positions ri and r2 at time t as^ 

A,(ri,r2,t;r',t') oc Vp(r',i') J d^ki J d^fca 

^ giki-(ri-r')-iwi(t-i')gik2-(r2-r')-iLJ2(i-t') ^2) 

In this paper, when the integration limits are omitted 
it is understood that the integrals are evaluated in the 
entire domain of the integration variables. 
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FIG. 1: Illustration of the nonlinear scattering of one photon 
into two by a scatterer at position r' at time t' . Photon 1 is 
observed in the position ri and photon 2 in the position r2, 
both at time t. 



In the above probability amplitude, we are coherently 
superimposing possibilities in which the energy and mo- 
mentum of photons 1 and 2 have all possible values, 
most of them whose sum is very different from the en- 
ergy and momentum of the incident photon, thus being 
possibilities that do not conserve energy or momentum. 
This treatment is closely related to the Feynman path- 
integrals formalism, and the conservation of energy and 
momentum in the process appears somewhat surprisingly 
as a consequence of the coherent sum of all possibilities. 
The point is that possibilities that do not conserve en- 
ergy, for instance, interfere destructively among them- 
selves. This will become clearer during our calculations. 

Let us now consider that instead of a single nonlinear 
scatterer we have a continuous set of nonlinear scatterers 
forming a nonlinear crystal, as depicted in Fig. [21 For 
simplicity, let us also consider that the linear response of 
the medium outside the crystal, represented by its refrac- 
tive index, is the same as inside, the nonlinear response 
being absent. In this situation, there is a fundamen- 
tal quantum indistinguishability of the position and the 
time in which the photon pair is created in the process 
of parametric downconversion. This condition forces us 
to coherently sum the probability amplitudes of gener- 
ating the photons in any part of the crystal and at any 
time to find the total probability amplitude of finding 
photons 1 and 2 in the positions ri and r2 at time t. 
Fig. [5] illustrates the probability amplitudes for generat- 
ing the photon pair in the spacetime points (ra,ia) and 
(rh,ff,). The twin-photon wave function '0(^)(ri,r2,i) is 
proportional to this total probability amplitude, and we 
can write 

V'(')(ri,r2,t) cx / dV / dt'A,{n,r2,t;r',t'), (2) 

where the integral in (fir' is in the crystal volume V and 
Ac(ri,r2,i;r',i') is given by Eq. IJ). 

Eqs. ([1]) and ([2]) are the essential equations of our 
treatment, and we summarize their physical interpreta- 
tion now. Eq. ^ indicates that the twin-photon wave 
function can be obtained from the coherent scattering 
of one pump photon into a photon pair in the nonlinear 
crystal. Since there is a fundamental indistinguishability 
of the position and the time in which the photon pair 
is created, the probability amplitudes for generation at 
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FIG. 2: Illustration of the nonlinear scattering of one pho- 
ton into two by a nonlinear crystal. Due to the fundamental 
indistinguishability of the position and time of the downcon- 
version process, all possibilities must be coherently added. In 
the figure the probability amplitudes of the conversion in the 
spacetime points {ra,ta) and {rb,tt) are represented. 



different times and positions must be coherently added, 
leading to the volume and time integrals in Eq. ([2]) . Fix- 
ing the time and the position of the twin photons cre- 
ation, the probability amplitude for this situation, given 
by Eq. ([ij, is computed taking into account the fact 
that, since the photons are created at precise time and 
position, they must have huge uncertainties on their en- 
ergies hui = hkic/rii and momenta ftki. All possibilities 
with different energies and momenta must be coherently 
added taking into account the corresponding phases accu- 
mulated during the propagation to the observation space- 
time points, leading to the integrals in dPki and d?'k2 in 
Eq. ID). The probability amplitude of the photon pair 
creation is also proportional to the pump photon wave 
function at the time and position of the conversion, as 
indicated in Eq. ([T]). 

In the next section we will find the twin-photon wave 
function evaluating the integrals in Eq. ([2]) using reason- 
able approximations. 



III. CALCULATING THE WAVE FUNCTION 
OF THE TWIN PHOTONS 

To calculate the wave function of the twin photons us- 
ing Eq. ([2|), we will first make some assumptions. If the 
incident pump photon is almost monochromatic, we can 
approximately write its wave function in position rep- 
resentation from the wave function in momentum rep- 
resentation as in the traditional quantum mechanics of 
massive particles 



^/'p(r',i') « / d^fcpVip(kp)e*''- 



(3) 



where we are writing the decomposition in terms of the 
wavevectors kp and angular frequencies ujp = kpc/up in- 
stead of the momenta Kkp and energies hup to simplify 

the notation. Decomposing kp 



•,2 ^ 



1p + ■\J f^p ~ Qp Qp 
being the component of kp in the xy plane, we note that 
Qp and ujp completely determine kp. So we can write the 
wave function in momentum representation in terms of 



Qp and Wp, changing the integrals in Eq. ([3]) from d^fcp 
to integrals in dqp and dujp: 



ipp{r,t)(x / dqp / dujp4:p{qp,ujp)e 



(4) 



This decomposition in terms of qp and uip is useful to 
treat cases in which the incident pump photon is in a 
beam mode that propagates nearly parallel to the z di- 
rection, in the paraxial regime. We will consider that 
^p(qp,a;p) has non- negligible values only for qp ^ fcp, 
such that the pump photon propagates nearly parallel 
to the z axis. Writing v' = p' + z'z, where p' is the 
component of r' in the xy plane, we have kp • r' = 



cip ■ p' + y kp ~ z'. Using a similar notation for the 

wavevectors of photons 1 and 2 and considering that the 
nonlinear crystal has dimensions in the x and y directions 
much larger than the width of the beam mode, such that 
the x' and y' integrals can be extended up to infinity in 
Eq. and a small dimension in the z direction, the 
wave function of the twin photons from Eq. ([2]) can be 
written as 

■0^^Hri,r2,t) (X J dqp J dujp J dqi J duJi J dq2 J duj; 
x4>p{qp,ujp) I dp'e'(i''-ii-'i^>''' / dt'e-'^""-"! 



^giki-ri-iwitgik2-r2-iwit 



u>2)t' 

(5) 



where the integral in z' gives a constant if we also have 
ki 3> qi and k2 ^ q2- In the appendix we discuss the 
relation between the omitted integral in z' and the phase 
matching condition for the efficient generation of photon 
pairs. The integral in p' results in a term proportional 
to the 2-dimensional delta function 6'^'^^ (qp — qi — q2). If 
the time light takes to propagate from the crystal to the 
observation points ri and r2 is greater than the duration 
of the incident photon pulse, the integral in t' can be 
extended until t' = oo, resulting in a term proportional 
to 5{ujp — cji — W2)- So the wave function of the twin 
photons can be written as 

V'^^'(ri,r2,i) oc j dqi J duJi J dq2 J duj2 

X i^p(qi +q2,o^i +a;2)e*^-""^-*"^*e*'^^-"-^-*"i*.(6) 

We can see that the same decomposition for a one- 
photon wave function used in Eq. ^ is used in Eq. 
([6]) for a two-photon wave function, such that the twin- 
photon wave function in the (q, w) representation is 

■!/i(^)(qi,wi,q2,W2) oc ?/'p(qi +q2,wi + ^2). (7) 

We can see in the above wave function that the sum 
of the energies of the twin photons {hioi + hw2) must 
be equal to an energy component of the incident photon 
{hujp), so energy is conserved in the process. Also, the 
sum of the momentum components of the twin photons in 
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the xy plane (ftqi + ftq2) must be equal to a momentum 
component of the incident photon (?iqp), so the transver- 
sal momentum is also conserved in the process.— Ac- 
cording to the appendix we also have hki^ -\-hk2z ~ fikpz, 
where kiz corresponds to the z component of the wavevec- 
tor of photon i, such that the total momentum is con- 
served in the process. But note that we did not impose 
conservation of energy or momentum in any part of the 
calculation. The energy conservation was a consequence 
of the coherent superposition of the probability ampli- 
tudes of the pair creation at any instant of time in the 
crystal, as the delta function S{LUp — ioi — UJ2) is a con- 
sequence of the t' integral in Eq. ([5]). In the same way, 
the conservation of momentum is a consequence of the 
coherent superposition of the probability amplitudes of 
the pair creation in any position inside the crystal, as 
the 2-dimensional delta function (5^^-'(qp — qi — q2) is a 
consequence of the p' integral in Eq. ([Sj. These conser- 
vation laws can be obtained simply from the interference 
of the probability amplitudes associated to indistinguish- 
able situations. 



IV. ENERGY AND MOMENTUM 
ENTANGLEMENT OF THE SYSTEM 

The twin-photon wave function of Eq. ([7]) cannot be 
written as a product of a wave function for photon 1 and 
a wave function for photon 2. This condition charac- 
terizes the entanglement of the system. If the incident 
photon has narrow energy (hujp) and momentum (Kkp) 
distributions, we see that the sum of the energies of the 
generated photons is equal to the energy of the incident 
photon, whereas the energy of each photon can assume 
a large range of values, characterizing a state which is 
highly entangled in energy. The sum of the momenta of 
the generated photons is equal to the momentum of the 
incident photon, whereas the momentum of each photon 
can also assume a large range of values, demonstrating 
that the quantum state of the twin photons is also highly 
entangled in momentum. 

An entangled state of the form of Eq. (O was used 
by Einstein, Podolsky and Rosen (EPR) in their famous 
paper on the completeness of quantum mechanics>i2 Let 
us consider the wave function for the x component of the 
photon momenta when the incident photon is a plane 
wave propagating in the z direction. In this case, we can 
write iI'^'^\pitP2) oc S{pi+p2), withpi = liQix- Of course, 
in a realistic case the delta function must be seen as an 
approximation for a very narrow distribution. On the 
other hand, if we consider the wave function in position 
representation at the exit face of the crystal (considered 
to be extremely thin), we have ip'-^^xi, X2) oc S{xi — 0:2), 
since the twin photons are born in the same position. 
Again, the delta function must be seen as an approxi- 
mation. Consider now that the generated photons reach 
two observers Alice and Bob who are placed in arbitrarily 
separated regions of space. If Bob measures the x mo- 



mentum component of his photon and obtain the value 
P, we can state with certainty that if Alice measures the 
X momentum component of her photon she will obtain 
the value —P. So EPR say that there is an dement of 
physical reality associated to the momentum of Alice's 
photon. On the other hand, if Bob measures the x posi- 
tion of his photon using lenses to project an image of the 
crystal in the observation plane and obtains the value X, 
we can state with certainty that if Alice measures the x 
position of her photon in a similar way she will obtain the 
value X, and EPR say that there is an element of phys- 
ical reality associated to the position of Alice's photon. 
However, if Bob is arbitrarily far from Alice, his measure- 
ment cannot affect her photon in any way, so EPR argue 
that Alice's photon should have elements of physical re- 
ality associated to both position and momentum. Since, 
in quantum mechanics, two observable quantities repre- 
sented by noncommuting operators can not have defi- 
nite and (in principle) predictable values simultaneously, 
EPR argue that quantum mechanics is not a complete 
theory, as "hidden variables" that are not considered by 
the theory would be essential to guarantee the simulta- 
neous reality for the position and momentum of Alice's 
photon in the gedanken experiment discussed above. 

Bohr's replyi^ to EPR argument was that the mea- 
surement of position or momentum made by Bob are mu- 
tually incompatible experiments. Bob's choice of which 
experiment he will perform on his photon determines dif- 
ferent types of prediction that he can make for experi- 
ments made by Alice. And there is no experiment that 
Alice can perform on her photon that could reveal what 
experiment was performed by Bob. Bohr argued that 
quantum mechanics is indeed a complete theory, since it 
can predict the probability of the experimental results of 
any combination of compatible experiments, which is the 
actual objective of the theory, such that there were no 
paradox. 

The discussion continued in philosophical grounds un- 
til the seminal work of BcU.^^ Bell showed that there are 
quantum entangled states with correlations between two 
parties that are stronger than what is allowed by "hidden 
variables" theories. The fact that the experiments con- 
firmed the quantum mechanical predictionsSS (although 
there are still some loopholes) definitely changed our way 
to see nature. We cannot have a consistent local and real- 
istic description of nature, since in an entangled state like 
the one from Eq. ([7]) we cannot attribute reality to each 
of two separated photons. We must consider either that 
Bob's measurement instantaneously change the state of 
Alice's photon, violating locality (but not permitting any 
instantaneous transmission of information between Bob 
and Alice) , or that the experimental results of the mea- 
surements performed by Alice do not depend only on the 
properties other photon (independently of Bob's photon) 
and on her measuring apparatus (independently of Bob's 
apparatus), violating any reasonable definition of reality 
for her photon properties. Accessible discussions of these 
fundamental aspects of quantum mechanics can be found 
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in the chapter 6 of Ref. d and in Ref. 

A physical implementation of the EPR state was per- 
formed by Howell and co-workers using the twin photons 
produced by parametric downconversion as described 
above. ^ The high degree of entanglement of the state of 
Eq. ([7]) makes the twin photons generated by paramet- 
ric downconversion to be an extremely valuable resource 
for testing this kind of fundamental aspects of quantum 
mechanics. The energy entanglement of the system leads 
to similar behaviours, since there is a complementarity 
between time and energy similar to the one between po- 
sition and momentum (although there is no time opera- 
tor in quantum mechanics). Moreover, entangled states 
have correlations among parties of the system that are 
stronger than what is allowed by classical physics, and 
some quantum information protocols take advantage of 
these quantum correlations to be more secure or efficient 
than the classical ones J The twin photons generated by 
parametric downconversion have also been extensively 
used in the implementation of these quantum informa- 
tion protocolsii 

V. CONCLUSIONS 

We derived an expression for the quantum state of 
the twin photons generated in the process of paramet- 
ric downconversion by a scattering procedure that coher- 
ently sum the probability amplitudes for the photon pair 
to be generated in any position inside the crystal and at 
any time, since these are fundamentally indistinguishable 
situations. The quantum state we find is equivalent to the 
one obtained with the traditional perturbative approach 
to calculate the Hamiltonian evolution of an electromag- 
netic field interacting with a nonlinear mediumi^^— Our 
treatment, however, is more intuitive, giving useful phys- 
ical insights to the problem. The state is highly entangled 
in both energy and momentum, making this system very 
useful for experimental tests of fundamental aspects of 
quantum mechanics and for the implementation of quan- 
tum information protocols. A more formal treatment of 
the work presented here can be found in Ref. 
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Appendix: Phase matching in parametric 
downconversion 

In this appendix we discuss the necessary phase match- 
ing for the efficient generation of photon pairs in para- 
metric downconversion. The phase matching condition 



can be obtained from the integ ral in z' omitted from Eq. 
([5]). Approximating kf — qf « ki — qf/{2ki), we have 



J~L/2 

(A.l) 

where the width of the crystal in the z direction is L. 
The efficiency of generation of photon pairs is propor- 
tional to /, so we must have kp = ki + k2 if we want 
a coUinear generation, in which the generated photons 
propagate close to the z axis, otherwise the oscillations 
of the first exponential above with z' will decrease the 
value of /. The wavenumber of each photon can be writ- 
ten as kiiuji) — ujiUiiuji)! where ni(y)i) is the refractive 
index for the photon i with angular frequency Wi. Let 
us assume that each of the photons is post-selected in 
a relatively narrow range of frequencies, such that the 
refractive index for each photon can be considered a con- 
stant rii. This can be obtained with the use of inter- 
ference filters in the photon detectors that are going to 
be used to detect the twin photons. Since according to 
Eq. ([7]) the sum of the frequencies of the generated pho- 
tons is equal to the frequency of the incident photon, 
the phase matching condition for coUinear generation can 
be written as rip = (ni -I- n2)/2. However, the inherent 
dispersion of the material makes the refractive index in- 
crease with frequency, such that in an isotropic medium 
we have rip > (n\ -f n^)/'^-, making it impossible to have 
an efficient generation of photon pairs. This problem is 
solved with the use of a nonlinear and birefringent ma- 
terial such that the refractive index also depends on the 
polarization of the photons, making Up — (ni -|-n2)/2 for 
photons with different polarization states. Type-I para- 
metric downconversion corresponds to the case in which 
both generated photons have polarization orthogonal to 
the polarization of the incident photon and in type-H 
parametric downconversion one of the generated photons 
has polarization orthogonal and the other one parallel to 
the incident photon. 

In cases where kp — ki — k2 = —a < 0, we can have 
an efficient generation of photon pairs that propagate in 
directions making angles with the z axis such that q\ /ki + 
q\/k2 = 2a and the second exponential in Eq. (jA.ip 
cancels the first, since Qp/kp is assumed to be negligible. 
For this reason the generated photons are usually emitted 
around cones. In any case we have kiz + k2z ~ kpz, where 
kiz corresponds to the z component the wavevector of 
photon i. 

Another important approximation that we make in our 
treatment is to consider that the integrals in c?qi and 
dq2 in Eq. ((5|) can be made without considering the 
dependence of I on qi and q2 in Eq. (jA.ip . This is 
a good approximation for crystals with small width L, 
since the accumulated phase in the second exponential 
in Eq. (jA.ip becomes negligible once the phase matching 
condition is established. 
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We could use spherical-wave modes to describe the propa- 
gation of photons 1 and 2 from (v' ,t') to {r\,t) and (r2,t) 
in Eq. ([l]). This would be more precise, but the plane- wave 
decomposition that we use, in which all wavevectors ap- 
pear with the same probability amplitude, leads to simpler 
calculations with essentially the same results. 
If the photon is not monochromatic the change of repre- 
sentation is not so simple, as discussed in Refs. [glliol. 
It is important to stress that the form ftk is not the only 
possibility for the momentum of a photon with wavevector 
k in a medium, as it is discussed in Refs. 15 16. There are 
many different ways of dividing the total momentum of 
an electromagnetic wave in a linear medium into electro- 
magnetic and material parts, all of them compatible with 
momentum conservation, and for each chosen division the 
photon momentum will have a different expression. The 
form ftk is the simplest one to treat the present case, that's 
why we've used it. 



